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ABSTRACT 

Aims. We aim to constrain the location, composition, and dynamical state of planetesimal populations and dust around the young, 
sun-like (G2 V) star HD 107146. 

Methods. We consider coronagraphic observations obtained with the Advanced Camera for Surveys (HST/ ACS) onboard the HST in 
broad V (/Ir » 0.6 /jm) and broad I (/I,. » 0.8 yum) filters, a resolved 1.3 mm map ohlained with the Combined Array for Research in 
Millimeter-Wave Astronomy (CARMA), Spitzer/IRS low resolution spectra in the range of 7.6yum to 37.0yum, and the spectral energy 
distribution (SED) of the object at wavelengths ranging from 3.5 //m to 3.1 mm. We complement these data with new coronagraphic 
high resolution observations of the debris disk using the Near Infrared Camera and Multi-Object Spectrometer (HST/NICMOS) aboard 
the Hubble Space Telescope (HST) in the Fl lOW filter (/l^ ~ 1.1 yum). The SED and images of the disk in scattered light as well as in 
thermal reemission are combined in our modeling using a parameterized model for the disk density distribution and optical properties 
of the dust. 

Results. A detailed analytical model of the debris disk around HD 107146 is presented that allows us to reproduce the almost entire 
set of spatially resolved and unresolved multi- wavelength observations. Considering the variety of complementary observational data, 
we are able to break the degeneracies produced by modeling SED data alone. We find the disk to be an extended ring with a peak 
surface density at 131 AU. Furthermore, we find evidence for an additional, inner disk probably composed of small grains released at 
the inner edge of the outer disk and moving inwards due to Poynting-Robertson drag. A birth ring scenario (i.e., a more or less broad 
ring of planetesimals creating the dust disk trough collisions) is found to be the most likely explanation of the ringlike shape of the 
disk. 
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1. Introduction 

Debris disks were first discovered via infrared excess emission 
associated with main sequence stars such as Vega from obser- 



& Chiang ' 2006; Wyatt 2008V However, studying the spectral 
energy distribution (SED) of the dust alone provides only weak, 
ambiguous constraints for its properties such as chemical com- 



vations with the Infrared Astronomical Satellite (IRAS; Aumann 



position, grain size, and spatial distribution ( Wolf & Hillenbrand 



|et al . 1984). The first spatially resolved image of such an object 
was obtained around/? Pictoris ( Smith & Terrile|1984| l. In the last 
few years the Spitzer Space Telescope has increased the number 
of known debris disk candidates identified by their thermal emis- 
sion above photospheric levels to several hundreds. As the most 
readily detectable signposts of other planetary systems, debris 
disks help us to improve our understanding of the formation and 
evolution of them as well as of our own solar system (Meyer] 
et al. |2007[ |Wyatt|2008 [ pCrivov|20 1 0| l . The dust detected in such 
disks is thought to be removed from those systems by the stel- 
lar radiation on short time scales compared to their ages. Thus, 
it must be transient, or more likely continuously replenished by 
ongoing collisions of bigger objects like planetesimals left over 
from the planet formation process. For isolated systems the cir- 
cumstellar dust's spatial distribution can be influenced by gravi- 
tational interaction with planets in addition to the influence of the 
star pominik & Decin|2003||Kenyon & Bromley|2004[|Strubbe| 



2003 ). For example, the location of the inner disk radius of the 
dust distribution is strongly degenerate with the lower limit of 
the grain size distribution. Direct measurements of the spatially 
resolved surface brightness by imaging mitigate many of these 
degeneracies. Images of disks seen near to face-on are particu- 
larly useful to directly determine the azimuthal and radial distri- 
bution of the dust and thus to trace the density structure. 

In this paper, we present new spatially resolved observations 
of the HD 107146 debris disk obtained with the Near Infrared 
Camera and Multi-Object Spectrometer (NICMOS) onboard the 
Hubble Space Telescope (HST). These data complement existing 
spatially resolved images at optical and millimeter wavelengths 
as well as spatially unresolved data of the star and the disk over a 
broad range of wavelengths (Table [T]i. An approach to combine 
all available data and so to derive a detailed model of the disk 
and to break degeneracies produced by modeling single datasets 
is presented. 
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Fig. 1. Scattered light images. The HST/ ACS images at 0.6 yum (left) and at 0.8 yum (center) and the //5r/NICM0S image at l.ljum 
(right). The images are oriented north-up and the HST/ ACS images are rebinned to the HST/NICMOS scale of 0.076 "/pix. The 
plotted coordinate system indicates the distance from the star in arcseconds along the two axes. Note the different flux scales in 
mJy/as^ and the fact that the stellar PSF residuals are heavily truncated in these scales. The object in the SW is a background 
galaxy. 



HD 107146 is a young, well studied, solar-type star (G2 V), 
28.5 ± 0.7 pc ( Ferryman et al.|1997 1 from the sun. Its luminos- 
ity and effective temperature are ^ 1.1 L0 ( Ardila et al.||2004) l 
and » 58 50 K, respectively. Its estimated age ranges from 80 to 
320 Myr ([Williams et al.|2004t IMoor et al.|2006[|Roccatagliata 



et al.||2009[ l. |Silverstone| ( |2000[ ) was the first to note the pres- 
ence of an IR excess based on IRAS measurements. The origin 
of this excess in a debris disk system was discussed by Metchev 
let al.| ( |2004"l l. HD 107146 has been well studied in the context 
of the Spitzer Legacy Science Program FEPS (Formation and 
Evolution of Planetary Systems; Me yer et al.|20 06) resulting in 
a broad variety of photometric measurements as well as spectra 
of the star and the disk (e.g. Hillenbrand et al.|2008) . The debris 
disk around HD 107146 is remarkable in many respects. With its 
large inner hole of ~ 60 AU and its extension of > 200 AU in 
radius, corresponding to 2" and 7", respectively ( [Ardila et al.| 



2004|l, it can be imaged by coronagraphy without masking large 



parts of the disk. Its particularly high scattered light surface 



brightness enables high signal-to-noise imaging ( [Ardila et al.[ 
[2004[ l. The disk is seen with an inclination of 25° + 5° from 
face-on ([Ardila et al.|2004ji. This means, that the radial surface 
brightness distribution can be observed without strong projec- 
tion effects which would occur in a close to edge-on case. With 
the inclination known, one can properly account for the remain- 
ing effects. In addition, the inclination of the disk enables us 
to constrain the dust properties from potential scattering phase 
asymmetries. 

We present our new observations in Sect. [2] The available 
data collected from the literature is summarised in Sect. [3] 
Details about our modeling, particularly the modeling software, 
the applied model, the order in which we include the available 
data and the information we extract in each step, respectively, 
can be found in Sect. |4] The results are presented and discussed 
in detail in Sect. l5] Conclusions are given in Sect. l6] 



2. New observations 

//577NICMOS coronagraphic observations of HD 107146 were 
obtained on 24 January 2005 as part of a larger circumstellar 



disk imaging surve}[J(i/5r GO/10177; PI: Schneider) designed 
to reveal and spatially map near-infrared starlight scattered by 
dust orbiting the central stars in both protoplanetary and more 
evolved debris systems (Sc hneider et al.|2005| l. Our HD 107146 
observations were conducted in our pre-planned survey mode, 
unchanged as a result of the discovery of a circumstellar scat- 
tered light disk around HD 107146 at 0.6 /im and 0.8 ;um wave- 
lengths with ACS coronagraphy ( [Ardila et al.|2004 1 shortly be- 
fore our already committed NICMOS observations. 

The HD 107146 field was observed with NICMOS camera 
2 in two contiguous HST visits at celestial orientations differ- 
ing by 29?9 in a single 54min HST target visibility interval. 
HD 107146 was placed behind the 0'.'3 radius coronagraphic hole 
using ACCUM mode target acquisition exposures (0.75 s) in 
the F187N filter at each visit. Three 224 s (STEP32/NSAMP14) 
multiaccum exposures were taken in the Fl lOW (A^s - 1 . lO/^m, 
FWHM = 0.59 /vm) filter Following the deep coronagraphic im- 
ages, short pairs of 8.53 s (SCAMRR/NSAMP21) FllOW ex- 
posures of the star were obtained with the star moved outside 
of the coronagraphic hole to directly enable flux density scaling 
with template point spread function (PSF) candidate stars simi- 
larly observed in GO/10177. In each visit, the individual instru- 
mentally calibrated count-rate images derived from the corona- 
graphic multiaccum exposures, were median combined and cor- 
rected for the well characterized 0.7 % geometrical distortion in 
NICMOS camera 2 (rectified to a two dimensional linear image 
scale of 75.8 mas pix"') yielding two images at different celes- 
tial orientations each with an effective integration time of 672 s. 
Additional technical details concerning the observational strat- 
egy and data reduction techniques employed for all GO/10177 
targets are presented in , Schneider et al., ( ,2005, ,2006^ . 

In our survey, in addition to our science targets, we simi- 
larly (but more deeply) observed five very bright stars of spectral 
types B through M (spanning the range of our targets) that were 
qualified through previous HST observations for use as coron- 
agraphic PSF subtraction templates^ including GJ517 (K5V; 
J=6.876) and GJ 355 (KOV; J=6.076). The set of five corona- 



' http://nicmosis.as.arizona.edu:8000/POSTERS/AASJAN_2006 
_GO10177.pdf 



see: http://www.stsci.edu/observing/phase2-public/10177.pro 



Steve Ertel et al.: Multi-wavelength modeling of the HD 107146 debris disk 



graphic PSF images created from these observations were co- 
aligned and flux-density scaled to HD 107146 as described in 
Schneider et al. (2005). Circumstellar scattered light excesses 
were seen at r > 2" in both fully reduced coronagraphic images 
of HD 107146 after separately subtracting each of the candidate 
PSF template images. These circumstellar structures, which con- 
form to the morphology of HD 107146 seen in scattered light at 
optical wavelengths as reported by Ardila et al. ( 2004| ) from ACS 
imagery, are absent in identically processed PSF minus PSF tem- 
plate subtractions (as shown, for example, in Fig. 1 of Schneider 
et al. ( 2006| l for two of the same PSF templates used here). 

The detailed sub-structures within the GJ517 and GJ355 
coronagraphic PSFs were most closely matched to those arising 
from the underlying stellar component in the HD 107146 im- 
ages, and when separately subtracted yielded the most artifact- 
free circumstellar disk images. This was not wholly unexpected 
since GJ 5 17 and GJ 355 are closer in |J - H| color to HD 107146 
than our other PSF candidate stars^ While GJ 355 is nearly the 
same brightness as HD 107146 (in the Fl lOW passband), GJ 517 
is significantly fainter (AJ = 1.0). This deficiency was partially 
mitigated by the longer and deeper integrations employed in our 
Fl low PSF template star observations (for GJ 5 17, 1 120 s total 
integration time) compared to our disk targets. 

The four PSF subtracted images of HD 107 146 using GJ 5 17 
and GJ 355 as PSF subtraction templates were median com- 
bined, after rotating all about the co-aligned position of the 
occulted star to a common celestial orientation, except in in- 
dividually masked and rejected regions significantly degraded 
by post-subtraction residual optical artifacts (primarily, insuffi- 
ciently suppressed diffraction spikes from the HST secondary 
mirror support). This resulting two-orientation re-combined im- 
age that reveals the HD 107146 circumstellar scattered light disk 
is shown in Fig. [T] along with the ACS images (Sect. [3]) that are 
used in the scattered-light analysis discussed subsequently in this 
paper. 

Because of the nature of the survey observations (with modal 
overheads for two-roll coronagraphy) allocating only a single 
spacecraft orbit per disk target, the NICMOS integration was 
rather shallow (22 minutes, combining all exposures at both field 
orientations) compared to ACS imaging in GO program 10330 
(H. Ford PI) at shorter wavelength (x 39min and 36min at 0.6 
and 0.8 /vm, respectively). The 1.1 yum two-roll re-combined PSF 
subtracted image of the HD 107146 circumstellar disk (Fig. [T] 
right panel) is shown in a bi-polar linear display stretch in sur- 
face brightness from -40juJy/as^ (hard black) to -1-120 //Jy/as^ 
(hard white) alongside the PSF subtracted images obtained with 
ACS as discussed by|Ardila et al. (2004). The circumstellar disk 
seen in the NICMOS image is instrumentally truncated in a small 
sector to the northeast at the apparent radius of the bright debris 
ring of a: 3'.' 5 due to the field-of-view limitations of NICMOS 
camera 2, but is clearly revealed at all other circumstellar az- 
imuth angles. Residual artifacts after PSF subtraction (empha- 
sized for visibility relative to the brightness of the debris ring in 
the image display for illustrative purposes with the stretch cho- 
sen) are manifest primarily in two forms. First, high spatial fre- 
quency radial streaks and, second, small amplitude zonal under 
(or over) subtractions at some radii (the mottled light and dark 
backgrounds at low spatial frequency), both which often arise in 
imperfect NICMOS coronagraphic PSF subtractions. Two ob- 



jects appear in the FOV, a galaxy « 6" from HD 107146 (to the 
SW) also seen in the ACS images, and a slightly diffused point- 
like image artifact (at Aa = +7'.'8, A6 - 0.0") persisting from 
the preceding target acquisition imaging. Neither of the PSF sub- 
traction artifacts significantly affect the visibility of the disk, but 
their appearance does somewhat degrade the fidelity of the im- 
age. 

3. Available data and previous work 

In addition to the i/577NICMOS observations we consider vari- 
ous already published observations. In this section these data as 
well as some more theoretical works including HD 107146 are 
described in chronological order. 

|Silverstone| ( |20 00) identified HD 107146 as a candidate for 
having a Vega-like disk based on IRAS excess at wavelengths of 
60 /im and 100/im. Metchev et al. (20041 performed 10 micron 



observations of HD 107146 using the Palomar 5 m Telescope. 
While the lOyum emission was found to be photospheric, they 
independently from |S ilverstone ( 2000) identified HD 107146 as 
a Vega-like candidate based on the already mentioned IRAS 
observations. They estimated a stellar age of 50 to 250 Myr 
and a minimum dust mass of 6.6 x 10"^ Mq assuming gray 
body grains. Different masses were also derived using blackbody 
grains (6.0 x lO^'^Mo) and ISM-like grains (6.0 x lO^^M©) il- 
lustrating the strong dependence of the derived mass on assump- 
tions of grain composition. 

Williams et al. ( 2004^ performed SCUBA observations of the 



HD 107146 disk at 450A'm and 850yL(m. They found the disk 
to be marginally resolved with an extent of 300 AU x 210 AU 
(10'.'5 X 7'.'4) at 450/im (beam size - 8'.'0) but to be unresolved 
at 850//m (beam size = 14'.'5). They also performed Keckll 
observations using the facility instrument LWS at 1 1 .7 yum and 
17.8 yum, but did not detect any extended structure. They also 
modeled the disk's SED using a single temperature modified 
black body and found clear evidence of a large inner hole (> 
31 AU in radius). Based on this modeling they estimated a dust 
mass of a 3 X 10"^ Mq. A detailed discussion of the age of the 
star (30 to 250 Myr) was given in WilUams et al . (2004). 

Ardila et al. (2004 2005| resolved the disk using HST/ ACS 
PSF subtracted coronagraphic observations in the F606W (broad 
V) and F814W (broad I) filters (Fig. [T). These data ai-e the ba- 
sis for our resolved modeling of the debris disk. They found the 
disk to be a broad ring with maximum vertical optical depth at 
130 AU from the star and a FWHM of 85 AU. They fitted elUp- 
tical isophotes to the disk images and found the disk inclined 
25° ± 5° from a pole-on orientation with the disk minor axis at a 



celestial position angle (east of north) angle of 58° ± 5°. Ardila 



^ All three of these stars were observed after a first in four year mid 
HST Cycle 13 focus/desorption compensation move of the HST sec- 
ondary mirror that additionally constrained PSF template suitability in 
our GO/10177 targets. 



et al.,(2004) also pointed out that the debris disk exhibits a sur- 
face brightness asymmetry in their ACS data, consistent with 
forward scattering of the star light by micron-sized dust grains 
given the inferred disk inclination. 

[Carpenter et al. ( 2005[ ) marginally resolved the disk at 3 mm 
using the Owens Valley Radio Observatory (OVRO) millimeter- 
wave array. They found the total flux density at 3 mm to be con- 
sistent with the previous work of Williams et al. (2004). The 
angular size of the disk in that image is consistent with that de- 
rived from all other data, but the structure seems to be much 
more elongated. Because of the low resolution and low signal to 
noise ratio (SNR) we do not use this image for detailed model- 
ing, but include its total flux density in our SED analysis. 

IKrivov et al. p008) compared results derived in their work 
from a collisional approach using a birth ring scenario to the 
SED of HD 107146. They found the SED to be reproduced by a 



Steve Ertel et al.: Multi-wavelength modeling of the HD 107146 debris disk 



Table 1, Photometric data of the HD 107146 system 



Wavelength 



Flux Uncertainty (Icr) Reference 



[mJy] 



[mJy] 



3.6 

4.5 

8.0 

10.3 

11.7 

13.0 

17.8 

24.0 

33.0 

60.0 

70.0 

100.0 

350.0 

450.0 

850.0 

1300.0 

3100.0 



1711.3 

1074.8 

384.4 

247.0 

175.0 

138.9 

85.0 

59.8 

86.7 

705.0 

669.1 

910.0 

319.0 

130.0 

20.0 

10.4 

1.42 



36.7 

24.7 

8.2 

22.0 

10.0 

8.5 

8.0 

2.5 

5.7 

56.0 

47.8 

155.0 

45.0 

12.0 

3.2 

1.5 

0.25 



References. (l)IHillenbrand et al.l([2008l; (2)IMetchev et al.|p004t; (3 ) 
[Wilhams et al. (20041; (4) Moor et al. (20061; (5) Colde r etal.H2009> ; 
(6) |Najita & WilliamS|^2005^ ; (7) |Carpenter et al.|^2005^ 



parent belt producing the dust at ^ 200 AU, too large compared 
to the available, spatially resolved data. This discrepancy was 
discussed in the context of roughness of Mie calculations due to 
the assumption of homogeneous spheres that might be a rough 
estimate of the shape of the dust grains. 

|Hillenbrand et al.| ( |2008j ) presented unresolved Spitzer mea- 
surements at 3.6, 4.5, 8.0, 13, 24, 33, 70, and 160 jum all of which 
are included in our SED analysis. The low resolution spectra in 
the range of 7.6 to 37.0 /im obtained with Spitzer/IRS and used in 
their work are available through the Spitzer archive. We obtained 
a low resolution spectrum over the whole range as a combination 
of the available SLl, LLl, LL2 and LL3 spectra and included it 
in our present study. 

[Corder et al.| ( |2009| l presented a resolved 1 .3 mm map, ob- 
tained using CARMA . In this image a ring-like structure is 
clearly seen. In comparison to the scattered light images of 
|Ardila et"ar] (|2004 ) they found a less extended disk with the ring- 
like structure concentrated at ~ 97 AU and two possible peaks 
in surface brightness at about the same distance from the star. 
A marginally resolved image at 350 /im is also presented. The 
1 .3 mm map is considered in our resolved modeling of the disk. 

Roccataghata et al. ( 2009| l computed some basic properties 



of the HD 107146 debris disk using its SED (see their paper 
for details). They found an inner radius of 10.2 AU, a minimum 
grain size of 8.6 //m and a dust mass of 3.3 x 10"^ M0 as best- 
fit values from their modeling. The inner radius as well as the 
minimum grain size reported by Roccatagliata et aL] ( |2009 1 are 
inconsistent with that expected from the resolved data and the 
blow out size of the system. These discrepancies have been dis- 
cussed in the context of grain segregation. They state that the 
modeled SED data predominantly trace the position of large (up 
to millimeter size) grains in contrast to the smallest grains (near 
the blow out size) traced by the HST scattered light images. 

The compiled photometric data is listed in Table [T] 



4. Analysis 

For our analysis of the debris disk we compare simulated data 
from a dust reemission modeling suite to the available spectral 
energy distribution and multi-wavelength imaging data. 



4.1. Model description 

The simulations of the dust reemission SED ar e done using the 
software tool DDS 



( |Wolf & Hillenbrand|2005l l. For simulations 
of resolved images our newly created tool debris is used. It 
allows to simulate images of arbitrary optically thin dust config- 
urations in thermal reemission as well as in scattered light using 
exact Mie theory. The dust properties derived from Mie scatter- 
ing are thereby computed with our software-tool mi ex ( |Wolf & 
Voshchinnikov|2004| included directly in debris. 

The stellar emission from HD 107146 is estimated by fitting 
a [Kurucz, (^1969) model to the 0.4 - 2.5 yum stellar SED. The fit 
gives Teff = 5924 K and \ogg - 4.48. A distance of the star 
from the sun of 28.5 pc is adopted for our simulations. The de- 
rived uncertainties on the model parameters represent only the 
formal errors from our fitting. They do not include uncertain- 
ties in the distance of the star. For our cuiTent simulations of the 
disk an analytical, rotationally symmetrical density distribution 
n(r) is employed (r. distance from the star). We limit on a sin- 
gle component model (no multi-ring structure) to minimize the 
degeneracies of our modeling. There is no obvious evidence for 
a multi-ring structure at this point of the work. A possible ex- 
tension of our model motivated by the shape of the Spitzer/\RS 



spectrum is discussed in Sect. 5.2.1 



4.1 .1 . Employed density distribution 

The simplest approach for the radial dust distribution is a power- 
law distribution n{r) oc r^" from an inner radius Ri„ to an outer 
radius /?out of the disk. As shown in Fig. |2] the observed radial 
brightness profiles derived from the HST/ACS images exhibit a 
decrease with increasing distance to the star, similar to a power- 
law distribution only in the outermost parts of the disk, outwards 
of a: 130 AU from the star. Further inward the surface brightness 
profiles flatten significantly. This is inconsistent with an increas- 
ing particle density. The profile derived from the F606W image 
even decreases with decreasing distance from the star inwards of 
~ 100 AU. Inwards of this distance the profile derived from the 
F606W image is more reliable then that derived from the F8 14W 
image (Sect. 4.3 1. In the following we motivate and describe the 
employed density distribution. 

We search for a simple analytical function which closely re- 
produces the disk density profile using a minimum of free pa- 
rameters. Such a distribution allows simultaneous fitting of all 
parameters and can easily be compared to any distribution found 
in other studies. The actual distribution will be better described 
by a large number of physical parameters, reflecting all physi- 
cal processes being responsible for it (e.g. the distribution of the 
planetesimals producing the dust, the effect of radiation pres- 
sure and Poynting-Robertson drag on the dust, mass and orbits 
of possible planets influencing the distribution by gravitational 
interaction). If a good fit is achieved, our empirical approxima- 
tion will at least be very similar to this distribution. 

As described above, the observed surface brightness profile 
can be described by an increase followed by a decrease with in- 
creasing distance from the star. In addition, the position of the 



http://www 1 .astrophysik.uni-kiel.de/dds/ 
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peak and the width of the distribution have to be parameterized. 
Not Hmiting ourselves on functions with (known) physical inter- 
pretation, we consider the following candidate density distribu- 
tions, that might be able to reproduce the described behaviour: 

- A product of two power-laws, one with positive and one with 
negative exponent: In such a distribution there is a strong 
correlation between the single parameters, making the fitting 
complicated and potentially resulting in very extreme, not 
physically interpretable values for all parameters. 

- A product of two exponential functions can be simplified to 
one single exponential function and does not allow to repro- 
duce the observed behaviour. 

- A product of a Gaussian distribution and a power-law: Due 
to the symmetrical Gaussian dominating the behaviour of the 
function, distributions with strong asymmetries around the 
peak distance can only be achieved by very extreme indexes 
of the power-law part. Moreover the parameters in this dis- 
tribution are heavily correlated making fitting complicated. 

- A product of a power-law and an exponential function (anal- 
ogous to Planck's law): As shown below, the number of pa- 
rameters and correlations between them can be reduced to a 
minimum for this distribution. 

Based on the above discussion we decide to reject the first 
three candidate distributions and to adopt the last one which is 
most convenient for our purpose. From an ad hoc point of view 
the following distribution can be employed: 

(0 V r < s y r < Tsub 

In this equation the quantity rjub is the sublimation radius of the 
dust, ro is a scale length stretching the distribution along r, s is 
a characteristic length shifting the distribution along r, and a^ 
and a2 are parameters describing the slopes of the dust distribu- 
tion in the inner and outer region, respectively. In this equation 
one can express the quantity ro through the peak radius rp of the 
distribution using the following condition: 



dn'(r) 
dr 



(r - r„) -0; r^ > s 



(2) 



Defining a new density distribution n(r) normalised to n(rp) - 1, 
one can eliminate the parameter a2'- 



n{r) = 




■ exp 



ffi 



I '■p-V 



V r < i V r< Tsub 
^ r^ s A r > Tsub 



(3) 



In the resulting density distribution there are only three free pa- 
rameters, s, Tp and a\, describing the inner increase and the outer 
decrease of the density as well as its width and peak position. It 
is important to note, that the number of free parameters is still 
the same as for the power-law distribution (a, Ri^, and Rom), de- 
spite an apparently much more complex shape of the distribu- 
tion. Furthermore, only two parameters, s and a\, are correlated 
(see discussion in Sect. 4.3 i, while even in the case of a simple 
power-law a, Rin, and Rom are strongly correlated because of the 
poor fit on the observed data. In addition, due to the particular 
shape of this distribution, there will be usually no need for an in- 
ner cut off Rin, which is essential in the power-law distribution. 

In debris disks we often observe an additional outer break 
in the surface brightness suggesting a break in the radial den- 
sity distribution (e.g. |Kalas et al.||2006[ ). Such a break can be 



realized in the same way as for the case of a power-law distri- 
bution employing an outer radius rout of the density distribution: 
n(f > '"out) = 0. While an outer radius of the disk is usually 
needed to limit the extension of the disk for a power-law distri- 
bution, that falls to negligible low values very slowly, the expo- 
nential decay reaches such values much earlier. So if there is no 
real, physical break observed in the surface brightness profile (as 
it is the case for HD 107146), the outermost radius will just be 
larger than the flux limited extent of the disk, and is only used to 
limit the extension of our simulated images. Thus, in this study 
we set rout to 500 AU. In contrast to the power-law distribution 
with the two cut-off radii our new distribution is continuous and 
continuously differentiable. Due to its analytical character it can 
easily be compared to any distribution found in other studies. 

Additionally, a constant opening angle of the disk of 10° is 
employed to define the vertical extension of our disk. That means 
the disk height is increasing linearly with distance from the star. 
As long as the angle is small enough, this assumption is of lower 
importance for our modeling. This simplification is chosen be- 
cause of the nearly face-on orientation of the disk and allows one 
to compute the surface density distribution E(p) from the radial 
density distribution n(r) very easily: 



E(p) oc n(r) ■ r 



(4) 



Here the quantity p is the distance from the star projected on the 
midplane of the disk. 

4.1.2. Employed dust grain properties 

The only available data that could inform about the mineral- 
ogy of the dust grains is the Spitzer spectrum. Since this spec- 
trum does not exhibit any significant features, we employ for the 
grain com position astronomical silicate with a bulk density o f 
2.7g/cm^ ( JDraine & Leel|l984[ [Weingartner & Draine||2001 1. 
The grain size distribution follows a power-law 



n{a) oc a 



(5) 



from a minimum grain size a^i^ to a maximum grain size amax- 
We assume the grain size distribution to be the same at all dis- 
tances, so we can compute the abundance «(r, a) of a certain 
grain size at a certain position in the disk through 



n(r, a) oc n(r) ■ n{a) 



(6) 



This assumption is necessary to limit the complexity of our 
model. It is not fully justified, since the effects of radiation pres- 
sure as well as Poynting-Robertson drag are stronger for smaller 
grains (larger fi ratio. Fig. l3]l which naturally leads to grain 
segregation. As the disk is expected to be collision dominated 
(Sect. 4.3 I, we assume that grain segregation has a minor effect 
on the distribution of grains of different size. This assumption 
will prove successful to reproduce most of the available data on 
HD 107146. 

In our grain size distribution Omin and y are free parameters 
but flmax is fixed to 1 mm, since the influence of bigger grains on 
the thermal reemission as well as on the scattered light can be 
neglected for significantly steep size distributions. Therefore the 
dust masses derived in this work represent the mass in particles 
smaller than 1 mm in radius only. 

4.2. General guideline 

In order to efficiently break the modeling degeneracies one has to 
understand the information that can be extracted from the obser- 
vations and how to combine the ancillary data. In the following 
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we present the general guideline used to reduce the number of 
free parameters to a minimum in each step of the modeling. At 
the end of the modeling study each of our free parameters de- 
scribed above has been fitted. The described approach leads to a 
complete set of model parameters for one unified model of the 
disk, consistent as far as possible with all included data. 

Radial dust distribution from high resolution images: Having 
high resolution scattered light images of a radially symmetrical, 
optically thin debris disk seen near to face-on, the surface num- 
ber density profile E(p) of the disk is correlated with the face-on 
surface brightness profile ip{p) as follows: 



Y(p) (X p' ■ ^(p) . 0(90°) 



(7) 



Because the scattering efficiency at a scattering angle of 90°, 
<I)(90°), of the dust is unknown at this point, one can not derive 
the absolute value of the surface density. For detailed fitting one 
has to make reasonable assumptions on certain dust properties 
(in our case a„ 



and y, see Sect. 4.3 



for the assumptions made 
in the present work). This is possible, because these assump- 
tions primarily affect the total brightness of the disk and thus the 
resulting disk mass but not the shape of the disk. The radial den- 
sity distribution n{r) can be found using Eq. |4] Due to the near 
to face-on orientation, the effect of the vertical structure of the 
dust distribution on the surface brightness distribution is negligi- 
ble. Because scattered light images trace particularly the small- 
est grains, the validity of the derived density distribution for all 
grain sizes has to be verified. Additionally, for coronagraphic 
scattered light images the inner parts of the disk are hidden by 
the coronagraph or heavily affected by the residuals of stellar 
PSF subtraction. Hence the density distribution in the inner re- 
gion remains uncertain. Both the validity of the distribution for 
larger grain sizes as well as in the inner region can be verified 
including spatially resolved observations at thermal reemission 
wavelengths, e.g. from mid-infrared to millimeter 

Dust properties from SED modeling: Once the spatial dis- 
tribution of the dust is constrained from modeling the scattered 
light imageajthe main degeneracy in SED fitting is broken. One 
can now derive properties of the dust grains, in our case the lower 
grain size and the exponent of the grain size distribution, from 
fitting the SED. One then has to go back to the scattered light 
images and confirm the results for the density distribution using 
these dust properties instead of the assumptions made before. 
Eventually, this process requires several iterations. 

A complete set of model parameters: With the final set of 
model parameters one can obtain (almost) independent mass es- 
timates from all images as well as from the SED. These masses 
should be consistent. In the case of a slightly inclined disk and 
scattered light data available one can further verify the fit and 
break potential remaining degeneracies quantifying how well the 
scattering asymmetry is reproduced by the model. 

4.3. Modeling the Hubble scattered light images 

For the two HST/ ACS images (Fig.[T]l |Ard"ilaetal.| ( p004l l found 
that the artifacts from the PSF subtraction dominate the scattered 
light within x 2" from the star. For a detailed discussion of the 
errors in these images see |Ardila et al.|p004| l. 

Fig. |2] shows the sky-plane deprojected (i.e., if seen face- 
on and assuming / = 25°) azimuthally medianed radial surface 
brightness profiles of the disk derived from the two HST/ ACS 



2. 0.01 




250 



Fig. 2. Azimuthally medianed (one pixel wide annular zones) 
surface brightness profiles derived from the HST/ ACS scattered 
light images after deprojection (for details see Sect. |4.3[ ). Note 
that the F606W radial brightness profile is reliable only for 
p > 60 AU and the F814W image is reUable only for p > 100 AU 
(dashed, vertical lines). An additional error of » 5% due to un- 
certainties in stellar PSF subtraction is not included in the plot- 
ted error bars. The stellar flux densities are 5.9 + 0.3 Jy and 
6.9 + 0.3 Jy at the F606W and F814W wavebands ( |Ardila et al.| 
j2004j). The black, solid lines represent the profiles derived from 
ouf final model of the disk (Sect. [4|. The albedo of our em- 
ployed astronomical silicate grains (Sect. |4TTT2| is 0.56 and 0.57 
at A - 0.6 /im and A = 0.8 /vm, respectively, for the smallest 
(2.5 /im) grains dominating the scattering. 



images. They are obtained following 'Weinberger et al. ( 1999| l. 
After deprojecting the images to face on, adopting the inclina- 
tion of the disk of 25° from [Ardila et al. (2004 1 and using flux 
conservation, the azimuthal median is derived in one pixel wide 
concentric radial zones. The error bar on each measurement is 
then computed following the equation 



1 



Y,(x-x,f 



(8) 



i=l 



^ Simplifying assumption: density distribution is the same for all dust 
species (size, composition; see Sect. 4.1.2 1 



where A^ is the number of pixels in a radial zone, x is the me- 
dian flux in a zone and xi is the flux in one pixel. This is anal- 
ogous to the standard error of the mean, but using the median 
in each zone. Deriving the azimuthal median does not account 
for the scattering asymmetries which are expected to be elimi- 
nated when medianing azimuthally. This assumption is less rel- 
evant because the same procedure has been applied deriving ra- 
dial profiles from the modeled images, so measured and modeled 
profiles can be compared directly. In the Fl lOW image large PSF 
subtraction residuals do not allow to derive a reliable profile. 

The disk surface brightness derived from the two HST/ACS 
images decreases with increasing radius for r > 130 AU. A simi- 
lar behaviour can also be seen qualitatively in the HST/NICMOS 
image. Within a radius of 130 AU of the star the two HST/ACS 
surface brightness profiles flatten significantly. Inwards of » 
100 AU the F814W image shows a behaviour that is completely 
different from that of the F606W image. While the F606W im- 
age shows a decrease in surface brightness with decreasing ra- 
dius, the surface brightness obtained from the F814W image 
exhibits a steep increase. An examination of the F814W data 
shows that the best global scaling of the PSF template results 
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in an under-subtraction of the stellar light at 2.0" < r < 3.5" 
(57 AU < r < 100 AU), particularly toward the SE. This be- 
comes obvious when comparing the F814W and F606W images 
in Fig. [T] Furthermore, no physical process is known to us that 
could produce such a different behaviour at 0.8 /vm compared to 
0.6 fim and 1 . 1 fim. Following the ACS Instrument Handbook, 
Version 8.0 (Boffi et al."2007). Sections 5.6.5, 6.2.9, and 9.3.2, 
and Version 10.0 ( Maybhate et al._2010j . Section 5.6.5, there is 
a known halo in the HRC images beyond 700 nm. Furthermore, 
the PSF in this filter is relatively large, so that less light is 
blocked by the coronagraph. This results in a coronagraphic PSF 
that is more sensitive to spot shifts and star-to-spot misalign- 
ments. Thus, we attribute the increase of the surface brightness 
profile derived from the F814W image with decreasing distance 
from the star at r < 100 AU to technical artifacts and residuals 
from imperfect PSF subtraction. 

Due to the lower amplitude of the PSF subtraction residuals 
compared to the other scattered light images, the F606W image 
is the best choice for starting our simulations. To reduce the in- 
fluence of PSF subtraction errors on our modeling, we mask the 
inner 2" (» 6QAU) in radius around the center of the image men- 
tioned by Ardila et al. (2004) and neglect it in the analysis. To 
obtain comparable results, the same region of our simulated im- 
ages is masked. An error map for the PSF subtracted image is 
obtained from under- and over-subtractions of the PSF in each 
pixel derived from different PSF reference stars observed. The 
fonowing equation is employed: 



1 



■ ( I hi; - low,- 1 ) + noise 



(9) 



In this equation / runs over all pixels and low, and hi, are values 
obtained from the under- and over-subtraction (minimum and 
maximum value in one pixel subtracting the different PSF ref- 
erences). Additionally, we add the Icr noise level derived from 
pixel to pixel noise in regions of the original image far from the 
star where no signal is detected. 

To find the radial dust distribution, the parameters describing 
the grain size distribution are fixed and only the free parameters 
of our particle density distribution (ai, r^, and s) are varied. The 
parameters of the grain size distribution (Omin, amax. and y) are 
fixed as follows: 



- The value of Omax is fixed to 1 mm (Sect. [4TL2| l. 

- The value of fl,„in is fixed at 1.0 jum. The /3 ratio, namely the 
ratio between radiation pressure force and gravitational force 
on a dust grain around HD 107146 is computed (Fig. [3]l. A 
blow out size (J3 - 0.5) of 0.5 /im can be found. In contrast, 
the dust is scattering significantly red with respect to the stel- 
lar spectrum in the range of 0.6 to 1 . 1 /vm (Fig. [T] |2]l, imply- 
ing the absence of small, Rayleigh scattering dust grains. For 
the adopted dust properties, only a minimal grain size larger 
than a: l.O/zm results in a red scattering disk in the observed 
range of wavelengths. 

- A value of y = 3.5 is chosen, being expected from an 
equilibrium coUision al cascade ( DohnanyTl|1969 l. Thebault 
& Augereau ( 2007) modeled collisional processes in de- 



bris disks and found collision time scales Tcoh much smaller 
than 10^ years for micron-sized particles, decreasing with in- 
creasing disk mass. The Poynting-Robertson timescale tpr 
for a particle with fl,nin = 1 yum to spiral onto the star from 
100 AU is calculated. A value of 1 .7x 10^ years can be found. 
Taking into account the high mass of the disk (Williams et al. 



2004 1, the expected collision time scale is about two orders 
of magnitude lower, e.g. 10^ years. Because of the decreas- 
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Fig. 3. /3 ratio (the ratio between radiation pressure force and 
gravitational force) for astronomical silicate with a bulk density 
of 2.7 g/cm^ around HD 107146. The dashed horizontal and ver- 
tical lines indicate /? = 0.5 and the corresponding blow out size 
of a = 0.5 /im, respectively. 



ing influence of Poynting-Robertson drag for larger parti- 
cles, we conclude that collisions dominate the evolution of 
the dust significantly (tcoii/tpr < 6.0 x lO"-'). 



For the adopted inclination of the disk of 25° from Ardila 
et al.'(2004) the adopted values of Oniin, flmax, and y produce an 
azimuthal scattering phase asymmetry in the surface brightness 
that is in good agreement with the observed one. 

To search for the best-fit model, we simulate scattered light 
images for a grid of model parameters (ai, r^, and s) in our an- 
alytical density distribution. After a first, coarse search a grid 
with smaller spacings centered on the best fit value of the coarse 
search is employed for our final fitting. The considered range 
of parameters and the employed spacings are listed in Table [2] 
Each simulated scattered light image is convolved with the corre- 
sponding PSF, derived with the Tiny Tim softwareFland is scaled 
to minimize the x^ derived from the observed and modeled im- 
age and the error map: 



X' 



-Z 



/ observation - model \ 






uncertamty 



-) 



(10) 



In this equation the sum runs over all pixels of the used images 
that have not been masked (Fig. pjl. The derived best-fit x^ for 
each model is used to evaluate the agreement between observed 
and simulated image. 

Fitting our density distribution based on the F606W image, 
we find best-fit values of ffj = 7.0+[|-^, rp = 120.3+J5-^ AU, and 
s - 25.0;^4Q AU, where the uncertainties are given as lo" confi- 
dence levels. These uncertainties represent only the formal errors 
from the fitting process. The distribution of the derived values 
of X^ is plotted in Fig.0 The plotted confidence levels are es- 
timated by computing the values oi x^ on a three dimensional 



http://www.stecf.org/instruments/TinyTim/tinytimweb 
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Fig. 4. Projection of the three dimensional ellipsoids of constant ;\f^ on the rp = const (left), s - const (center), and ai - const (right) 
planes from fitting the F606W image. The contour levels show Icr, 2cr, and 3cr levels of confidence, respectively. The cross indicates 
the position of the best-fit values. The plots illustrate very well the correlation between ai and s and the very weak correlations 
between these parameters and rp, respectively. For a detailed description of the parameters shown see Sect. 4.1.1 



Table 2. Parameter space for modeling the HST/ACS F606W 
image (upper, middle) and the SED (lower). 



Parameter Grid 



Range 



Spacing 



r, [AU] 
s [AU] 



coarse 



0.0- 
112.0- 
-10.0- 



14.0 
129.0 
-60.0 



0.5 
1.0 
1.0 



r, [AU] 
5 [AU] 


fine 


4.5 - 10.0 

118.2-122.7 

10.0-50.0 


0.1 
0.3 
1.0 


Qmin [fJ^m] 

y 


— 


0.5 - 20 
2.0-5.0 


0.1 
0.1 



grid in the parameter space and by estimating the probability p 
for each model j using the equation 



P^^^v[-{x]-xL)] 



(11) 



Using the derived density distribution, one is able to repro- 
duce the behaviour of the radial surface brightness profiles of the 
F606W image and outside ^ 100 AU of the F814W image, not 
affected by the PSF subtraction errors described earlier in this 
section (Fig.l2]). 

4.4. Deriving a^m and y from tine SED 

In the previous section, we derived the radial dust density dis- 
tribution from the scattered light images. Under the conditions 
described in Sect. |4.3| the results of this approach are indepen- 
dent from the grain properties such as lower grain size and size 
distribution. Employing this density distribution and assuming 
that it is the same for all grain sizes we can now simulate the 
SED of the debris disk to derive best-fit values for the parame- 
ters flniin and y of our grain size distribution. We start with the 
values assumed in Sect. I4.3l and do not limit ourselves on a cer- 
tain range of the parameters. However, since we do not expect 



values too far from what we assumed in Sect. 4.3 we effectively 
search the range of fl,nin and y listed in Table |2| For each fit the 
corresponding disk mass is derived and compared with the mass 
derived from the F606W image using the same model parame- 
ters (dust composition and radial distribution). 



We do not find any solution giving total disk masses from 
SED and scattered light that differ by less than a factor of 1.4, 
where the masses derived from scattered light are in any case 
larger than that derived from the SED. Taking into account that 
these mass estimates have been derived from observations which 
potentially trace different fractions of the dust distribution (e.g., 
grains of different size), and the uncertainties in flux scaling that 
arise from photometric calibration and stellar PSF subtraction, 
the resulting mass range is rather narrow and the derived masses 
are considered to be consistent. However, we are not able to 
reach our initial goal to derive one mass value from all data sets 
to make it one final model parameter as it would have been from 
fully simultaneous modeling. This discrepancy is discussed in 
Sect. 5.2.2 For this reason we decide to fit the SED employ- 
ing the dust density distribution derived from the scattered light 
data, but without the constraint that the derived masses should be 
identical. Instead, we derive independent disk masses from each 
data set for our final, unified model with the radial dust distri- 
bution derived from the F606W image and the dust composition 
derived from the SED of the system. 

Furthermore, we are not able to reproduce the Spitzer spec- 
trum and the rest of the SED data simultaneously (Fig.|6]l. Thus 
the Spitzer spectrum is only considered including the synthetic 



photometry derived from it at 13 and 33 /im (Hillenbrand et al. 
|2008[ ). The discrepancies with the shape of the spectrum are dis- 
cussed in Sect. l5.2.TI 

Keeping in mind the above caveats we continue fitting the 
data with our approach described in Sect. |4.2| Employing the 
radial dust density distribution derived in Sect. |4.3| and the ap- 
proach to fit the SED described above, the best-fit values for the 
parameters describing the dust properties are amin = 2.5 /zm and 
y = 3.6 (Fig. l6]l. The SED can be only poorly reproduced us- 
ing flmin < 2.0 fim or Amin > 3.0 yum and any value of y as well 
as using y < 3.5 or 7 > 3.7 and any value of a^in- These val- 
ues are considered as limits of the fitted values. Adopting these 
results for the grain properties instead of the assumptions made 
in Sect. |4.3| we were able to confirm the results on the density 
distribution derived there. The value of y = 3.6 + 0.1 is consis- 
tent with that expected from an equilibrium coUisional cascade 
(r = 3.5; |Dohnanyi[T969| . 
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Fig. 5. Observed (top), modeled (center), and model subtracted images (bottom) in the F606W (left), the F814W (center), and the 
FllOW (right) HST filters using the final model derived in this paper. The images have been rotated, so that the major axis of 
the disk is oriented along the x-axis. Only the considered parts of the images are shown, while regions of large residuals from 
PSF subtraction or contaminated by other objects on the line of sight are masked. In the F606W and F814W images the ling-like 
structure of the disk is well reproduced by the model, and only the residuals from PSF subtraction of the central star remain visible. 
In the residual Fl lOW image parts of the disk remain visible in the residual image suggesting that scattering at this wavelength is 
underestimated in the model. Note the different flux scales in different columns. The model image in the right column is displayed 
in a different way than the data and model subtracted data to better visualize the disk. 



4.5. Verifying the model using the CARMA 1.3mm map. 

The scattered light images do not allow to determine the dust dis- 
tribution within the inner 2" (57 AU) from the star. Furthermore, 
the spatial dust distribution seen in scattered light is dominated 
by the distribution of the smallest grains and can put only weak 



constraints on the distribution of grains with radii larger than 
~ lOOjum. To verify the spatial dust distribution derived from 
the scattered light images one can use the CARMA 1.3 mm map 
(Fig. [Til. However, in this image most of the disk is detected with 
an average surface brightness of only x 3cr while there are two 
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peaks reaching ^ 5(t. The value of cr has been derived in regions 
of the image where no signal from the disk is expected. Due to 
its low SNR and the coarse angular resolution relative to the HST 
images we do not use the 1.3 mm map to fit model parameters. 
However we verify that our model is consistent with the CARMA 
observations. 

Therefore, we simulate images at 1.3 mm using the model 
derived above. We simulate observations of the best-fit model 
image using the actual u-v spacings from the CARMA observa- 
tions. These model u-v data points were subtracted from the ob- 
servations, and the residual image was reconstructed using the 
same image parameters as in [Corder et al. ( 2009[ l. To correctly 
scale the modeled data to the observed source flux, we minimize 
the remaining flux in the reconstructed, model subtracted image. 

From the reconstructed, model subtracted image (Fig. |7]l 
we find that our model reproduces the global behaviour of the 
CARMA map, even though better results for slightly lower peak 
radii and/or a shallower decrease of the density distribution in the 
inner regions might be found. While a small amount of residual 
flux remains visible in the inner region of the disk, in the outer 
regions the disk flux is over-subtracted. However, since nega- 
tive flux values with a ringlike shape around the disk are already 
visible in the original image and the flux remaining in the inner 
region of the disk hardly exceeds 3cr, we rate these deficien- 
cies as insignificant. We conclude that our density distribution 
is consistent with the CARMA data, keeping in mind the large 
uncertainties. Furthermore, the two peaks visible in this image 
are subtracted very well by our radially symmetrical image and 
peak residuals are around 3cr. This implies that these structures 
may be noise, rather than real structures in the disk induced by a 
massive planet orbiting the star within the inner edge of the disk 
as suggested by Corder et al. ( 2009| l. 

At this part of the work we can already conclude that the 
small particles seen in scattered light and the larger particles 
traced by the CARMA data are cospatial within the errors of the 
available data. This is important, since it supports the assump- 
tion of our model, that particles of all sizes are distributed in the 
same way. 



>, 100 - 




10 100 

wavelength [/^m] 



1000 



Fig. 6. Simulated SED from our single component best-fit 
model. The solid line represents the modeled SED, while the 
dashed line represents the stellar photospheric flux. Photometric 
measurements in the plotted wavelength range with correspond- 
ing error bars and the considered Spitzer spectrum are plotted, 
too. The discrepancies, in particular with the Spitzer spectrum 
are discussed in Sect. l5.2.TI 



CARMA data has been described in Sect. |4.5| The slight overes- 
timation of the disk extent compared to these data may explain 
the larger total disk mass derived from them compared to the 
mass derived from the SED. 

The derived disk masses are listed in Table [3] The scatter 
in the estimates based on the different data sets gives an idea 
of the self-consistency of our modeling approach with the data. 
Changing any of the model parameters within their confidence 
levels incurs correlated changes in some of the other parameters 
to maintain the quality of the fit (i.e., see Fig.|4]i, and affects the 
estimated disk masses by a factor of two to three. These changes 
in the disk mass are then correlated among the particular sets of 
data (scattered light images, SED and CARMA map). 



4.6. Fitting fhe final model to the available data and deriving 
disk masses 

In the above sections we employed the radial density distribu- 
tion of the dust described in Sect. 4. 1.1| (Eq. l3]l as well as the 
grain size distribution described there (Eq. l5]l to derive a full set 
of model parameters (ai, rp, s, a^in, and y). We can now de- 
rive a total mass of the disk from each data set. This is done 
for the SED, the F606W, F814W, and FllOW images by scaling 
the modeled data and subtracting them from the observations to 
minimize the residuals. We mask the inner 60 AU in the F606W 
and FllOW images, where the errors dominate significantly the 
signal. While the scattered light from the disk is subtracted very 
well in the F606W image, large parts of the disk remain visi- 
ble after subtracting our best-fit model from the FllOW image 
(Fig. |5]l. This seems to be due to a different distribution of the 
scattered light flux, closer to the star in this image, which re- 
sults in a scaling factor for fitting the model to the image that is 
too small to reproduce the total flux. However, due to the large 
uncertainties in this image mentioned earlier, we do not further 
consider this fit. In the F814W image the inner 100 AU affected 
by the large errors from PSF subtraction (Sect. |43] l and two ad- 
ditional image artifacts with large negative values that affect sig- 
nificantly the fit are masked. The remaining disk structure can 
be reproduced very well (Fig. |2] B). The approach used for the 



5. Discussion 

5.1. Results 

We used one model (with the total disk mass derived from the 
different data sets varying by a factor of » 1.5) to reproduce 
the surface brightness in the F606W images and the outer parts 
of the F814W image being not affected by the PSF subtraction 
residuals discussed in Sect.|43](Fig.|5]|2]), the CARMA 1.3 mm 
map (Fig.|7]i and the SED of the HD 107146 debris disk (Fig.|6]l. 
We find the disk to be a broad ring with a peak of the 
density distribution n(r) at 120.3^^5 AU, which corresponds to 
131.4+°g AU for the surface density distribution 2(r) (Fig. [sjl. 
The radial density distribution and surface density distribution 
transform into each other following Eq. |4] Using the derived 
confidence levels of the parameters, we find the FWHM of the 
density distribution to be 86;^p AU which corresponds to an 
FWHM of the surface density distribution of 91^}' AU. Th is is in 



good agreement with the values derived by Ardila et al. ( 2004 
2005 ). We find a lower dust grain size of 2.5^Qj/zm, inconsis 
tent with the blow-out size of the system of 0.5 jum. The expo- 
nent of the grain size distribution of y = -3.6;^q[ is consistent 
with the value expected from an equilibrium coUisional cascade 
( |Dohnanyi| 1 969| l . Table [3] summarizes the derived model param- 
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Fig. 7. Observed CARMA 1.3 mm map (left), modeled image (center), and model subtracted map (right). The contour levels have 
increments of Icr. White, dashed levels start at -Icr and indicate negative values, while black, solid levels start at Icr and indicate 
positive values, where Icr - 0.35 mJy/beam. A bold line is used to mark the maximum contour in each image. The FWHM of the 
synthesized beam (3'.'2 x 2'.'7) is shown in the lower-left corner. The images have been rotated, so that the major axis, as derived 
from the scattered light images, is oriented along the x-axis. 




Table 3. Derived Model Parameters 



100 150 200 250 

r [AU], /} [AU] 



350 



Fig. 8. Volume density distribution n(r) over the radial distance 
r from the star and surface density distribution E(p) over the dis- 
tance p from the star projected on the disk midplane for our best- 
fit model (parameters see Table |3]l. The analytical functions for 
n(r) and 2(p) transform into each other following Eq. l4] Both 
distributions are given in arbitrary units and scaled to a peak 
height of one for each distribution so that they can easily be 



compared. The peak particle density is 8.03 x 10 m , while 
the peak surface density amounts to 2.64 x 10^ m"^. 



eters. For the SED as well as for each resolved image we derive 
a total dust mass, presented in the table, too. 

The available data is consistent with all dust in the system 
having the same radial distribution. This can be expected from a 
collision dominated disk ( Wyatt 2005 ). We find from integrating 
the synthetic stellar spectrum and the modeled SED of the disk a 
fractional luminosity L(just/istai - 1.07x10"^ which supports the 
scenario of a heavily collision dominated disk. However, there 
are some data available that are reproduced by our model only in 
an unsatisfactory way, in particular the shape of the SpitzerflRS 
spectrum. These discrepancies are discussed in the following 
section. 



Parameter 


Derived Value 


Uncertainty* 


tti 


7.0 


6.2 ... 7.6 


rp [AU] 


120.3 


119.8... 120.7 


s [AU] 


25 


21.0... 33.0 


flmm [j"m] 


2.5 


2.0... 3.0 


"max [^m] 


1000.0 


(fixed) 


7 


3.6 


3.5 ... 3.7 


A/f606W [Mq] 


6.2 X 10-' 


- 


A^F814W [Mq] 


6.5 x 10-' 


- 


A^FllOW [Mq] 


8.5 x 10-' 


- 


McARMA [Mq] 


6.7 X 10-' 


- 


MsED [Mq] 


4.4 X 10-' 


- 



Notes. The uncertainties are given as Icr confidence levels. For details 
see Sects. 1431 |4!4l and|4.6| 



5.2. Remaining discrepancies between our model and 
available data 

We are able to reproduce most of the available data with only 
one global model of the disk. However, there remain some dis- 
crepancies that are discussed in this section: 

- At wavelengths A < 40 //m the modeled SED poorly re- 
produces the observed one sampled by the Spitzer spectrum 
(Fig.|6|. 

- We find a minimum grain size that is by a factor of 5 larger 
than the expected blow out size. 

- Residual flux remains visible at peak level of ^ 3cr in the 
inner region of the model-subtracted CARMA 1.3 mm map, 
while a slight over-subtraction is visible in the outer region 
of the disk. This implies a slightly lower extent of the disk 
at this wavelength and results in too large an estimate of the 
total disk mass from this image. 

- Significant parts of the disk remain visible in the Fl lOW im- 
age after subtracting our best-fit model, while we already 
derived a mass that is by a factor of x; 1.3 larger than that 
derived from the other scattered light images. 
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- We were not able to derive a unique value for the mass from 
all data sets. In contrast, the derived masses differ by a factor 
of up to 1.5. 

While the model allows to reproduce the overall structure of 
the disk very well, important insights may be gained from the 
remaining discrepancies. They may be ascribed to a number of 
deficiencies of our model that become only relevant being faced 
with such a large number of high quality complementary data: 

- A discontinuity of the radial density distribution in the inner 
regions of the disk hidden by the coronagraph in the scattered 
light observations (e.g., an additional, inner disk component) 
would be able to significantly alter the shape of the SED at 
the short wavelength edge of the excess. 

- The size distribution of the smallest grains may be only 
poorly described by a single power-law that describes the 
distribution of the larger grains very well. More detailed 
grain size distributions can be derived by detailed dynamical 



and collisional modeling of debris disks (e.g. Thebault et al. 



2003 



et al. 



Krivov et al.|2006||Thebault & Augere au 2007 ; Kiivov 
2008|l. The results show a wavy distribution with an 



overdensity of small grains (a; twice the expected blow out 
size). 
- A more complex chemical composition of the dust than the 
used astronomical silicate may be present. This has been 
studied for debris disks with significant emission features 
in the mid-infrared (e.g. Beichman et al. 2005; Esse et al.| 
[2007 2008 i C hen et al.||2008| l. A large variety of chemical 
dust compositions has been found, such as crystalline py- 
roxene, crystalline oUvine, amorphous olivine, and water ice 
inclusions. For disks without such features the effect on the 
whole SED has been studied (Miiller et al.|2010| Augereau et 
al. 20 11 in preparation). Introducing porous grains instead of 
the compact, spherical grains used in this study also changes 



the optical properties of the dust (e.g., Voshchinnikov et al. 
|2006| . 

- The derived density distribution as well as the surface bright- 
ness profiles derived from the scattered light data may be 
eiTor-prone due to limitations from the PSF subtraction. This 
might result in eiTors for the derived radial density distribu- 
tion of the dust and thus affect all derived model parameters. 

- The spatial distribution in our model is the same for all grain 
sizes. In contrast, dynamical modeling of debris disks shows 
that several influences can result in a more complex depen- 
dence of the dust distribution on grain size (e.g. Wyatt 2005 , 
[2006). The influences of radiation pressure and Poynting- 
Robertson drag, that increase with decreasing grain size (in- 
creasing j6-ratio), may effect the distribution of the smallest 
grains which results in a segregation of grains of different 
sizes. Thus, a narrow ring of planetesimals can produce a 
broad disk through collisions (birth ring scenario; |Strubbe &[ 
rChiang„2006l [Krivov et al.|[^06l |Krivov|[20T0| . However, 
since we find the disk to be colUsion dominated, we do not 
expect grain segregation to be a dominant process for this 
disk (|Wyatt|[2005 ). Even the presence of a massive planet 
does not result in strong radial grain segregation, although 
different spatial structures can be expected for different grain 
size ( |Wyatt|2006| l. 

In particular the effect of the first two points might have a 
strong effect on the shape of the SED at the short wavelength 
edge of the excess. Hence, we discuss possible extensions of 
our model in the following to account for these deficiencies. 
Exploring the chemical composition of the dust, in particular the 



shape of the dust grains, requires significant additional compu- 
tational effort and is beyond the scope of the present paper. 

5.2.1 . Unsatisfactory fit to tine Sp/fzer/IRS spectrum 

The most obvious and significant discrepancy of our model with 
the data is the different shape of the modeled SED in the range 
of the Sp/fzer/IRS spectrum. This might be explained in different 
ways. A more detailed grain size distribution than a power-law 
distribution as discussed above might result in a better fit to these 
data. Our model is able to account for such changes in the over 
all grain size distribution only by changing the power-law expo- 
nent and the lower grain size of the applied distribution, which 
might be insufficient. Since small (warm) grains dominate the 
emission at the short wavelength end of the excess, such an over- 
abundance of small grains is expected to result in an increase of 
the simulated flux in the range of the Spitzer spectrum. In such 
a scenario, one would expect that the additional excess flux seen 
in the Spitzer spectrum originates from the same region of the 
disk like the long wavelength excess. 

An alternative explanation for the discrepancy between our 
model and the Spitzer spectrum is an additional, inner disk 
present in the system, such as discussed in the case of the 
eEridani debris disk (Backman et al.|2009) . Such an additional, 
warm disk component would add excess flux at mid-infrared 
wavelengths and thus would be able to explain the missing flux 
from our model in this wavelength regime. In this case, the long 
wavelength and short wavelength excess should originate from 
different regions in the disk. 

None of the available, resolved data allow us to distinguish 
between the two scenarios. The HST data do not allow to con- 
strain the dust distribution in the innermost regions of the disk 
due to the coronagraphic nature of these observations. The com- 
parison of the model with the CARMA data gives a hint for addi- 
tional emission in the inner region of the disk, but the sensitivity 
of these data is too low, so that the deviation between model and 
data is insignificant. Furthermore, the disk is unresolved by the 
SpitzerlMiPS observations. 

Our approach to distinguish between the two scenarios, is 
to explore what one would have to add to reproduce the Spitzer 
spectrum using 

(a) An additional grain size component to mimic a more com- 
plex grain size distribution. Here the dust is spatially dis- 
tributed in the same way as the other dust (Sect. |4.3| l. The 
chemical composition of the dust is the same as before (as- 
tronomical silicate). We use only one distinct grain size here, 
so this is the only additional, free parameter. 

(b) An additional, inner disk component. We use a disk model 
that is described by a power-law radial density distribu- 
tion n{a) oc a"^ , ranging from an inner radius R'^'^ to an 
outer radius R'^^- The grain size distribution is described by 
n{a) oc a^^ ' , ranging from a lower grain size a'^'^^ to a fixed 
upper grain size of 1 mm (the same parameterization as the 
grain size distribution in our original model). We use a fixed 
opening angle of 10°, so that the vertical disk structure is 
described in the same way as in Sect. |4.3| The same chemi- 
cal composition of the dust as before (astronomical silicate) 
is employed. Since there is a total of 5 free parameters here 
and none of them can be fixed by simple physical assump- 
tions, it is clear that this fit will be very degenerate. However, 
we will be able to evaluate, whether this scenario is able to 
reproduce the shape of the Spitzer spectrum at all. 
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First, we change the parameters of our (outer disk) model so 
that the simulated flux of the system is smaller than the mea- 
sured flux in the Spitzer spectrum, but is still consistent with 
the long wavelength SED observed. Since our radial density dis- 
tribution derived from the scattered light data is still valid, we 
can only change the parameters of the grain size distribution. 
We find an increase of the lower grain size a^i^ from 2.5 /^m to 
3.5 fim without changing the exponent of the distribution, but in- 
creasing the total disk mass derived from fitting the SED data 
from 4.4 x 10"^ M© to 4.8 x 10"^ M© well suited to reach this 
goal. We then subtract the modeled SED of the system from the 
Spitzer spectrum and fit the remaining excess flux with the two 
additional disk components described above. The exact choice 
of flniin is rather arbitrary, since we are not able to distinguish 
between flux from the original disk model (with changed dust 
parameters) and the hypothetic additional component. Thus, it 
is not possible to derive any formal uncertainties on these pa- 
rameters. However, if we find an additional component that is 
able to reproduce the Spitzer/IRS spectrum, then in this ap- 
proach the value of a^in - 3.5 //m can be considered as a lower 
limit on the lower dust grain size of the original dust compo- 
nent. Smaller values would violate the constraint, that the excess 
from this component shall be smaller than the excess measured 
by Spitzer/IRS. On the other hand, almost any value larger than 
«min - 3.5 //m will be possible, since the missing flux at longer 
wavelengths can then be reproduced in a wide range by changing 
the parameters of the second, additional component. 

We find from our approach (a), that it is not possible to fit the 
shape of the Spitzer spectrum with an additional grain size com- 
ponent (with the radial density distribution derived in Sect. |4.3| l. 
Furthermore, this would add significant flux at longer wave- 
lengths. In any case the dust is too cold (too far from the star) 
in this approach. 

In contrast, with our approach (b) we are able to improve 
the fit on the Spitzer spectrum significantly without lowering the 
quality of the fit on the other SED data. The explored parame- 
ter space for our search for an additional, inner disk component 
and our best-fit results are shown in Table |4] Fig. |9] shows the 
new best-fit SED for this approach. Our best-fit gives a disk that 
is extended from the inner rim of the outer disk to very close 
to the star. The best-fit inner radius is identical with the small- 
est value of this parameter explored. This suggests that the disk 
is extended even to a lower inner radius. We find a nearly con- 
stant radial surface density distribution (exponent of radial den- 
sity distribution a'^'^^ = 0.8 which corresponds to an exponent 
of the radial surface density distribution of -0.2) to be the best- 
fit. The lower dust grain size (3.3 jum) is very close to the lower 
grain size of the outer disk in this approach (3.5 ;um). However, 
we find a very large exponent of the grain size distribution (best- 
fit: -y^'^^ = 10), i.e. a very narrow size distribution around the 
lowest grain size. The disk mass for this inner component then 
results to 3.6 x 10"" M©. From this model we predict total scat- 
tered light fluxes from this disk component at 0.6/im and 0.8/im 
of 0.037 mJy and 0.052 mJy. The total reemitted flux from this 
component is plotted in Fig. [9] as a function of the wavelength. 

The spatial distribution of the dust as well as the steep grain 
size distribution and the lower grain size that is consistent with 
that in the outer disk suggest, that the inner disk component in 
this scenario is produced by collisions close to the inner rim 
of the outer disk and then dragged inwards towards the star by 
Poynting-Robertson drag. This scenario has proven successful to 
explain the origin of the dust in the habitable zone of e Eridani 
( Reidemeister et al.|201 1 1. For a transport dominated disk with- 
out perturbation by a planet, a constant radial surface density dis- 



Table 4. Parameter space explored and fitting results for model- 
ing the remaining Spitzer/IRS flux by an additional, inner disk 
component. 



Parameter 



Range Spacing Best-fit Uncertainty' 



RT [AU] 


0.2-250.0 


temp- 


0.2 


0.2. 


.0.6 


Rf^ [AU] 


2.0-250.0 


temp" 


42.2 


6.0. 


.250 


a'"" 


-5.0-5.0 


0.1 


0.8 


-2.1 


.. 1.1 


<L L"m] 


0.5 - 100 


log-^ 


3.3 


0.8. 


.4.4 


T,add 


2.0-10.0 


0.1 


10.0 


3.6.. 


10.0 


M"" [Mq] 


free 


cont'' 


3.6 X 10-" 







Notes. ( 1 ) Confidence levels (Icr); (2) equally distributed in temperature 
of the dust species with the largest gradient of the radial temperature 
distribution (spacing = 2K); (3) equally distributed in log(a) (100 size 
bins); (4) continuous - The fitting of this parameter has not been done 
on a grid, but a continuous range of possible values has been used. 



1000 ^ 



^ 100 ^ 




10 100 

A [microns] 



1000 



Fig. 9. SED fit with the modified dust composition in the outer 
disk and an additional, inner disk component (Sect. 5.2.1 1. The 



two dotted curves indicate the emission from the two dust com- 
ponents only. 



tribution has been found there from dynamical modeling. Such 
a disk should consist predominantly of small grains. 

However, with the large uncertainties on our best-fit values 
it is not possible to draw strong conclusions from the resulting 
model. An additional, inner ring of planetesimals producing new 
dust through collisions (power-law exponent of the grain size 
distribution of 3.6) would also be possible within the \cr uncer- 
tainties of our fit. It is important to note that these uncertainties 
are only the formal uncertainties of our fit on the remaining flux 
in the Spitzer spectrum after subtraction of the SED from a suit- 
able model of the outer disk. They do not include the degrees 
of freedom introduced by the selection of the parameters of the 
outer disk which are expected to dominate the real uncertainties 
on the fitted model parameters. For example, a further increase 
of the lower dust grain size in the outer disk would be possible 
and would make a less steep grain size distribution more proba- 
ble. 

5.2.2. Different disk masses derived from different data sets 

While we were able to reproduce most of the available data using 
one model of the disk, we were not able to derive a unique disk 
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mass. The derived masses from the different data sets differ by a 
factor of up to 1 .5 (ignoring the HST/NICMOS data which show 
a larger deviation). In fully simultaneous modeling of all data, 
this mass would have been a parameter of the model that would 
have been to be fitted on the whole data set. This would have 
resulted in unsatisfactory fits on nearly all data sets, because of 
inadequate flux scaling of the different simulated data sets. The 
discrepancy is rather low, in particular keeping in mind the pho- 
tometric uncertainties in the single measurements, uncertainties 
from imperfect PSF subtraction, and the fact that we compare 
values derived from scattered light data and thermal reemission 
data. However, a global model of the disk should be able to re- 
produce all data without the need of different flux scaling done 
by deriving different masses from different data sets. 

Again, a more complex grain size distribution than the ex- 
pected one is one possible explanation. A different slope of the 
grain size distribution for the smallest grains present in the sys- 
tem would result in a different surface brightness at scattered 
light wavelengths, since these observations are very sensitive to 
the smallest grains. The slope of the grain size distribution in our 
single power-law model is mostly constrained by the shape of the 
SED at longer wavelengths (A > 100 ;um, which is most sensitive 
to larger grains. A derived disk mass from scattered light data 
that is too large (as it is the case for our results) would mean a too 
low surface brightness of our modeled data at these wavelengths. 
An underabundance of the smallest grains in our model would be 
a plausible explanation. However, we found in Sect. 



al. 2011, in preparation). As discussed before, adding smaller 
grains would in general also increase the scattering efficiency 
of the dust, increasing the modeled surface brightness and thus 
decreasing the mass derived from the scattered light images to 
a level more consistent with the value derived from fitting the 
SED. A different chemical dust composition or shape might also 
increase the value of ^6 of the dust and thus increase the expected 
blow-out size to a level more consistent with the findings of our 
modeling. Again, a different structure of the grains might have a 
very similar effect. 



5.3. Is the HD 107146 debris disk unusually massive for its 
age? 

Compared to other stars harboring debris disks, HD 107146 is 
in an intermediate age of « 100 Myr to 150 Myr. It is harbor- 
ing a very massive debris disk. Such high dust masses are much 
more usual for young debris disks around A type stars. Spatially 
resolved debris disks with comparable dust masses are for ex- 



amplejePic ( [Zuckerman & BeckUn|[T993|), HD 99803B (Wyatt 
et al."2003'), HR4796 dSheret et al.||2004| , and HR8799 (|Su 
et al. 2009). All these stars are A type stars or earlier. Of these 
stars HD99803B and HR 8799 might have a comparable age to 
HD 107146, while the two other ones are significantly younger. 
The dust mass in debris disks is expected to decline with age 



5.2.1 that 



adding an additional grain size component is not able to solve the 
discrepancy with the Spitzer spectrum. Furthermore, the small- 
est grain size is well constrained by the short wavelength end of 
the excess for the adopted chemical dust composition. Even with 
an additional inner disk component (Sect. [5^2T i there is a strong 
lower limit for the smallest grains. Thus, a different chemical 
dust composition, e.g., with larger albedo to increase the scatter- 
ing efficiency of the dust, is the most plausible explanation. 

An alternative to a different chemical composition of the 
dust could be a different structure of the grains, i.e. the pres- 
ence of fluffy, porous grains (e.g., |Voshchinnikov et"ar 2006 



trough well explored depletion mechanisms (e.g. |Silverstone 
200^|Lohne et al.|2008[|Roccatag liata et al.|2009|l and towards 



later spectral types due to lower initial disk mass. In contrast, 
HD 107146 is an intermediate age G type star harboring a disk 
with a mass comparable to the mass of disks around young 
A-type stars. This suggests that the HD 107146 disk is special 
among the known, spatially resolved debris disks. Its high age 
and high dust mass may imply that the system is undergoing 
an analog to the Late Heavy Bombardment (delayed stirring; 
Dominik & Decin 2003|l in our ow n solar system as also sug- 
gested for HR8799 lSu"et al.|2009| l. 



IVoshchinnikov & H"erm ing 2008|. This might alter the absorp- ^ ^ y^g oriqin of the disk shape 



tion efficiency and scattering efficiency and the wavelength de- 
pendence of these quantities independently. 

5.2.3. A large lower dust grain size 

The lower grain size from dynamical and collisional models of 
debris disks is in any case consistent with the expected blow out 
size, while there is an overdensity of particles with x twice the 
blow out size (e.g. Thebault et al. 2003 ). In contrast, we find a 
lower grain size of ~ 5 times the expected blow out size for our 
original model and even larger for the model with the additional 
inner disk component. Because small grains are in general hotter 
than large ones, they are traced particularly by the short wave- 
length end of the IR excess SED. While the lower grain size «„,;„ 
is usually heavily degenerate with the spatial dust distribution in 
pure SED modeling, this is well constrained in the present work 
by the scattered light images. As the smallest grains are the most 
efficiently scattering ones, our radial density distribution derived 
from scattered light images is expected to be particularly reli- 
able for these grains. Thus, grain segregation is not expected to 
have a significant effect on the lower grain size in our model- 
ing approach. To lower the value of flmi,, in our derived solution 
one can solely change the chemical composition or shape of the 
dust grains to get in general cooler dust. The inclusion of wa- 
ter ice has been successful to solve this problem (Augereau et 



The shape of the HD 107146 debris disk is similar to that 



|et al."200r), although with 6.2 x lO'^M© ( jSheret et al.||2004| i 
Fomalhaut's debris disk has a lower mass by about one order of 



of Fqmalha ut (HD 216956), an A3V star o f 200 Myr dHabing 

' ■■— ^' ■ — — «-• aLlP^^ 

by about one order of 
magnitude. Both stars have comparable ages. The similar shape 
of the two disks and the fact that the shape of the Fomalhaut disk 
seems to be induced by the planet Fomalhaut b orbiting the Star 
( Kalas et al.| |2008) raises the question if a planet is embedded 
also in the disk of HD 107146. [Kenyon & Bromieyl ( [2004) found 
from numerical simulations, that bright ring-Uke structures in 
debris disks are a by-product of planet formation while Thebault 
|& Augereau| (2007) simulated the collisional evolution of de 
bris disks without perturbing planets and without initial ling- 
like structure (e.g. a birth ring scenario) and found no evidence 
for the formation of ring-like structures or inner gaps. However, 
while the inner edge of Fomalhaut's debris disk represents a 
sharp break in particle density, the inner edge of the HD 107146 
debris disk seems to be a shallow decrease. Furthermore, the 
Fomalhaut disk is off-centered from the stellar position which 
is interpreted as a result of the planet-disk interaction. There is 
no obvious evidence in the data that the HD 107146 disk is off- 
centered in a similar way. Alternatively, a birth ring scenario of 
a collision dominated disk (as it is expected for the HD 107146 
debris disk) is also able to explain inner gaps in debris disks 
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without the need a massive interior planet ( |Krivov|2010| l. Here 
the planetesimals do not necessarily have to be distributed in a 
narrow ring, but the radial broadness of the planetesimal belt has 
to be similar to or smaller than that of the observed dust disk. 
Furthermore, the peak surface density of the planetesimal belt is 
expected to be close to the peak position of the surface density 
of the dust disk {^ 130 AU for HD 107146). 

Further //STcoronagraphic observations have been approved 
for execution to be conducted as part of the Cycle 1 8 program 
GO 1222^ (G. Schneider, PI), where 8 orbits of STIS coron- 
agraphy in the broadband optical are devoted to probe the inner 
regions of the HD 107146 disk as close as r = 0.3" (8.6 AU) 
and to look for resonant structures in the dust distribution. These 
data may, through higher quality imaging, better constrain pos- 
sible disk sub structures induced by planets as well as on the 
reliability of our derived density distribution, particularly in the 
inner regions of the disk. 



6. Conclusions 

A detailed model of the debris disk around HD 107146 has been 
created and fitted to all available data using the approach de- 
scribed in the paper A broad range of resolved data from optical 
scattered light to millimeter data as well as the SED has been 
reproduced by a single radial density profile. All parameters of 
our model have been fixed without remaining degeneracies us- 
ing only one scattered light image and the well-sampled SED 
of the disk in the particular case (Sect. |4]). We found the disk to 
be a broad ring of dust with smooth inner and outer boundaries, 
and with a peak position from the star and FWHM of the surface 
density of a; 131 AU and « 9 1 AU, respectively. We concluded 
from our modeling results that the disk is heavily collision dom- 
inated. Discrepancies of our model to the data are discussed. 
From these discrepancies we find strong evidence for an addi- 
tional inner disk component, possibly near the habitable zone 
(the stellocentric distance where liquid water can exist on plan- 
etary surfaces, i.e., few AU from a solar-type star) of the star 
and for a different chemical composition or shape of the dust 
grains. The additional inner disk component is predominantly 
composed of small grains and the results from our modeling are 
consistent with these grains being released at the inner edge of 
the disk and dragged inwards due to Poynting-Robertson drag. 
Only detailed, simultaneous modeling of the large amount of 
available high quality, complementary data was able to reveal 
these signposts. We did not find evidence of an orbiting planet 
from the available data and found that a birth ring scenario is 
likely responsible for the ringlike shape of the disk. 
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